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Abstract: Experiment NA49 at the Cern SPS uses a large acceptance detector for a 
systematic study of particle yields and correlations in nucleus-nucleus, nucleon-nucleus 
and nucleon-nucleon collisions. Preliminary results for Pb+Pb collisions at 40, 80 and 
158 A-GeV beam energy are shown and compared to measurements at lower and higher 
energies. 



1 Introduction 

The primary purpose of the heavy ion programme at the CERN SPS is the search for 
evidence of a transient deconfmed state of strongly interacting matter during the early 
stage of nucleus-nucleus collisions [|I|. The transition from a dilute state of individual 
hadrons to a phase of quasi-free quarks and gluons, the quark gluon plasma (QGP), was 
first suggested using qualitative arguments @ and later confirmed by quantum chromody- 
namics (QCD) on the lattice provided the energy density reaches sufficiently high values 
in an extended volume. 




Figure 1: Schematic layout of the NA49 experiment at the Cern SPS showing beam 
detectors, superconducting dipole magnets, time projection chambers (VTPC,MTPC), 
time-of-flight arrays (TOF) and calorimeters (RCAL,VCAL). A thin solid target T is 
used for A+A collisions (a), which is surrounded by a detector of slow protons (CD) for 
p+A collsions (c). A liquid H 2 target is employed for p+p collisions (b) 

The results from the heavy-ion program at CERN in fact indicated that deconfinement 
may set in within the SPS energy range 0. Within most model scenarios the data imply 
that the initial energy density exceeds the critical value. Originally proposed signatures Q 
of the QGP were observed in Pb+Pb collisions at the top SPS energy, i.e. J/^f suppression, 
strangeness enhancement, and possibly thermal photons and dileptons. The significance 
of these signals as QGP signatures has come under renewed scrutiny. Moreover, there 
is no observational evidence for a sharp phase transition from QGP to hadrons such as 
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phase coexistence 0, critical fluctuations |5|, || or more speculative effects like DCC |7|] 
or parity violation ||. 

The NA49 experiment is performing an energy scan from 20-158 A-GeV in an effort 
to strengthen the evidence for the onset of deconfinement by searching for anomalies in 
the energy dependence of experimental observables. Data at 40, 80 and 158 A-GeV have 
so far been recorded and analysed. A brief description of the apparatus (section 2) is 
followed by preliminary results for 7T, K, A and A yields, ttti correlations, event-to-event 
charge fluctuations, anisotropic flow (sections 3-6), and conclusions (section 7). 

2 The NA49 detector 

The NA49 experiment |J was designed for the investigation of hadron production in the 
most violent Pb+Pb interactions at the Cern SPS. The main features (Fig.ffl) are large 
acceptance precision tracking and particle identification using time projection chambers 
(TPCs). The first two are located inside superconducting dipole magnets which provide 
the particle trajectory bending necessary for momentum determination. Charged particles 
in the forward hemisphere of the reaction are identified from the measurement of their 
energy loss dE/dx in the TPC gas (accuracy 3 -4 %). At central rapidity the identification 
is further improved by measurement of the time-of-flight (resolution 60 ps) to arrays of 
scintillation counter tiles (TOF-T) and strips (TOF-G). Strange particles are detected 
via decay topology and invariant mass measurement. The forward calorimeter VCAL 
measures the energy of the projectile spectators from which one can deduce the impact 
parameter in A+A collisions. 




Figure 2: Left: rapidity distribution of 7r~ in central Pb+Pb collisons at 40 (squares), 
80 (triangles) and 158 (dots) A-GeV. Open symbols show values reflected at y*=0 (NA49 
preliminary). Right: total pion multiplicity (n) produced per wounded (participant) 
nucleon versus the Fermi energy variable F w s% 2 ^ for p+p reactions (open symbols) and 
central nucleus-nucleus collisions (full symbols). 
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3 Yields of tt, K, A and A 



Raw K + and K~ yields were extracted from fits of the distributions of dE/dx and TOF 
(where available) in narrow bins of rapidity y and transverse momentum py. The resulting 
spectra were then corrected for geometrical acceptance, losses due to in-flight decays and 
reconstruction efficiency. The latter was determined from embedding simulated tracks 
into real events and amounted to ~ 95%. Spectra of tt~ mesons were derived from the 
acceptance corrected negatively charged particle yields in pr and y (assuming the tt mass) 
by subtracting the estimated contribution of K~, p and the contamination from secondary 
hadron decays. The ratio tt + /it~ was determined in the region where both dE/dx and 
TOF are available (0.91, 0.94 and 0.97 at 40, 80 and 158 A-GeV) and was assumed to be 
y independent. 
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Figure 3: Rapidity distribution of K~ (top row) and K + (bottom row) from dE/dx 
(circles) and combined dE/dx and TOF (squares) analysis. Open symbols show values 
reflected at y*=0. (NA49 preliminary) 

The resulting rapidity distributions of ir~ are plotted in Fig.|2], left. The integrated 
yields are 312±15, 445±22 and 610±30 at 40, 80 and 158 A-GeV respectively. Pions are 
the dominant produced particle species and thus their number provides a measure of the 
entropy in a statistical model description of the reaction. The yield of pions (estimated 
here as (tt) = 1.5-((7r~) + (7r + ))) divided by the number of wounded nucleons (participants) 
Njy is shown versus the Fermi energy variable F = (-y/i — 2m N ) 3 / 4 / ^snn 1 ^ 4 ~ s^w i n 
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FigfJ, right. While p+p data show a linear rise throughout there is a change for A+A 
collisions (illustrated more clearly in the inset) in the SPS energy range. Below one finds 
a regime of slight suppression, above a region of enhancement with a steeper linear rise 
than for p+p reactions. This steepening has been interpreted as indicating the activation 
of a large number of partonic degrees of freedom at the onset of deconfinement JTH . 




Figure 4: Midrapidity ratio of K~/7r - (left) and K + / 7r + (right) as function of energy from 
NA49 (squares, preliminary) compared to measurements at lower and higher energies. 
Predictions of the RQMD (dotted) and HSD [12] (full curve) models are shown. 

Kaons contain about 75% of the s, s quarks in the produced hadrons at SPS energies 
and thus their number indicates the total strangeness content of the final state. The 
rapidity distributions are displayed in Fig.^| and integrate with small extrapolation to 
total yields of 18±1, 29±2, 50±5 for K" and 56±3, 79±5, 95±9 for K+ at 40, 80, 160 
A-GeV respectively. Yields of most particles, of course, increase with energy. Changes 
in the composition of the produced system are better characterised by particle ratios. 
Measurements of K~ /ir~ and K + /tt + at midrapidity from NA49 are plotted versus energy 
in Fig.f| and compared to results at lower and higher energies. 

A continuous increase is seen for Kr /n~ . For K + /tt + one finds a steeper rise followed 
by a maximum at the lower end of the SPS energy region and a gradual decrease. Within 
a reaction scenario based on nucleon-nucleon collisions these features might be attributed 
to thresholds and the decrease of the baryon density with increasing energy. The ratio 
K~ /tt~ exhibits the threshold of the KK production mechanism. The lower mass threshold 
of associate KY production leads to a steeper rise in K + /n + and the rapidly falling 
baryon density may result in a compensation of the declining contribution from the KY 
by the growing contribution from the KK production mechanism. However, the continuous 
increase seen for nucleon-nucleon collisions (e.g. Fig.|7], right) does not really support such 
an interpretation. 

The most abundantly produced hyperons are A and A for which Fig.|5| shows the 
rapidity distributions. These show a broad shape for A reflecting the associate production 
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Figure 5: Rapidity distribution of A (left) and A (right). Open symbols show values 
reflected at y*=0. A yields at 80 and 158 A-GeV are displaced vertically by 10 respectively 
20 units for clarity (NA49 preliminary) 



mechanism and the partially stopped participant nucleon distribution. In contrast the 
distribution is of narrower Gaussian type for A which are most likely produced as members 
of hyperon-antihyperon (YY) pairs. The ratios (A)/(7r) and (A)/(tt) of 4tc integrated 
yields are displayed as a function of energy in Fig.|| One observes a steep threshold rise 
for the (A) / (tt) ratio followed by a decline which can be mainly attributed to the rapidly 
decreasing net baryon density. The ratio (A) / (n) exhibits a continuous rise due to the 
high mass YY threshold. 
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Figure 6: 4n yield ratios (A) / (tt) (left) and (A)/(tt) (right) versus energy from NA49 
(squares, preliminary) and lower energy AGS experiments. The dotted line shows a 



prediction from the extended hadron gas model [13 



6 



Microscopic dynamical and statistical models have been used extensively to describe 
particle yields in a wide variety of reactions. The first class is often based on string 

followed 



excitation, fusion and hadronisation (e.g. HSD |TJ, RQMD (TJ, UrQMD [|Tj 
by reinteractions of the formed hadrons (RQMD, UrQMD). As seen from Fig.^ HSD does 
not reproduce the energy dependence of the midrapidity K + /tt + ratio. On the other 
hand, RQMD correctly predicts the trend, but somewhat overpredicts the ratio in the 
SPS energy range. 

A comparison with the Atc ratio (K + )/(tt + ) is presented in Fig.|7] left. Both UrQMD 
and RQMD get the steep threshold rise. UrQMD values are much too low in the plateau 
due to an overprediction of pions. RQMD does not follow the drop in the SPS range 
which is indicated by the measurements. 
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Figure 7: Left: Air ratio (K + )/ (tt + ) versus energy compared to predictions of the RQMD 
|TT| (dashed), UrQMD |L4| (dotted) and extended statistical Jl3| (solid curve) models. 
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Right: the strangeness content measure E s versus the Fermi energy variable F^ 
compared to the prediction of the statistical model of the early stage [lCj (curves). 

Since antihyperon production rates are small and isospin symmetry approximately 
holds ((K + ) « (K )) nearly half of the s quarks in the produced hadrons are contained in 
K + mesons. Moreover, strangeness conservation requires (s) = (s). Thus (K + ) measures 
to a good approximation one quarter of all s and s quarks in the final state hadrons. The 
energy dependence of the (K + )/ (tt + ) ratio (see Fig.^ left) therefore indicates a maximum 
in the fraction of strangeness carrying particles in the lower SPS energy range. 

Statistical models have been surprisingly successful in describing ratios of particle 
yields in many types of reactions over a wide energy range. Since the widths of ra- 
pidity distributions depend on particle mass in the SPS energy range and below, the 
model should preferentially be compared to 4tt yields. Fits of this model to NA49 data 
[I5fl have indicated that while there is relative hadro-chemical equilibrium in the strange 
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particle sector there seems to be an overall undersaturation with respect to non-strange 
particle yields. The statistical model as such makes no prediction concerning the energy 
dependence of particle production. However, it has recently been supplemented jTJ] by 
a parameterisation of the energy dependence of its two main parameters, baryochemical 
potential fis and temperature T. Predictions are shown in Fig.[7| left for (K + ) / and 
Fig.|6] right for (A)/ (tc). The trend of the data is well reproduced by this extended statis- 
tical model calculation. In detail the decrease in the SPS energy range of both (K + ) / (tt + ) 
and (A) / (tt) and even more so of (S) / (71*) (not shown) is not well described. 

Predictions have also been published for a statistical model which explicitly assumes a 
deconfmed phase in the early stage above a certain threshold energy |H| . In this model the 
strangeness to entropy ratio is assumed to be established initially and to persist through 
the hadronisation stage. A measure of this quantity is the ratio E s = ((A) + (K + K))/(tt) 
evaluated from the final hadron multiplicities which is plotted in Fig.^ right. After a rise 
corresponding to the purely hadronic reaction in the model, one observes a saturation 
at the a level consistent with the strangeness to entropy ratio expected for an initially 
deconfmed system. 
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Figure 8: Gaussian radius parameters R S ide, Rout, Riong fitted to the ir~7r~ correlation 
function evaluated in the longitudinally comoving frame plotted versus the average trans- 
verse momentum of the pair. Rapidity range y* < y < y*+0.5 (NA49 preliminary). 
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4 tttt correlations 



Correlations of pions with near equal momenta pi, P2 provide information on the size 
and internal dynamics of the fireball source at freezout The analysis was performed 
in the longitudinally comoving reference frame, decomposing the momentum difference 
Q = px - p 2 into long, side, out components. For 7r + 7r~ pairs the correlation peak at 
small Q is predominantly due to the Coulomb attraction and can be well reproduced by 
a Coulomb wave calculation jl7| using the measured effective source size. Correlations 
of 7r~7r~ are described by a product of parameterised Coulomb repulsion [l7| and the 
quantum statistics enhancement, fitted with a Gaussian parameterisation. The resulting 
radius parameters R S ide, Rout, Riong are plotted in Fig.| versus the average transverse 
momentum of the pair. No striking energy dependence is observed. Within expanding 
source models the decrease with Kt is a manifestation of the strong longitudinal {Ri on g) 
and radial {R S ide, Rout) flow at SPS energies. Moreover, there is little change in the 
lifetime To ~ Ri ong ^T / 'Mt or the emission duration At 2 « (Ro^-R-Lfe) / 'Pt °f the source. 
The small value of At does not indicate a soft point of the matter equation of state of 
the kind discussed in ref. [El near the onset of deconfinement. 



i 6r 1 

Q : MOGeV • 80 GeV Al60GeV 




resonance gas 



■ 



frozen QGP fluctuations 



Ay 




0.1 0.2 0.3 0.4 0.5 0.6 



< N ch>acc I < N ch>tot 



Figure 9: Measures of event-to-event charge fluctuations in central Pb+Pb colli- 
sions. Left: D versus the size of the rapidity window Ay. Right: $ g versus the ratio 
(N c h) I {N c h)tot of the multiplicity in the acceptance window and the total multiplicity in 
the events; the curve shows the prediction for independent particle emission plus global 
charge conservation (NA49 preliminary). 
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5 Event-to-event charge fluctuations 



Recently it was proposed that event-to-event fluctuations of the charge ratio R=N + /N_ 
or the net charge Q=N + -N_ may be sensitive to deconfmement in the early stage of 
nucleus-nucleus collisions [|1^, IjJ. The smaller charge quanta in a partonic phase are 



expected to result in a reduction of such fluctuations. 

The fluctuations of the charge ratio were investigated via the measure D [pD| which is 
corrected for the residual net charge in the considered rapidity interval Ay as well as for 
global charge conservation. The preliminary NA49 results are shown in Fig.|9] left and are 
found to be close to the expectation for independent particle emission plus global charge 
conservation (D « 4) and do not change significantly with energy. No evidence is seen 
for the reduction predicted for a resonance gas nor for the large decrease expected for 
a QGP phase. It is, of course, not clear whether reduced fluctuations in the QGP will 
persist through the hadronisation, rescattering and resonance decay stages. 

The quantity $ q was proposed in ref. ^Tf for studying fluctuations of the net charge. 
It is independent of the number of superimposed particle sources, has value zero for 
independent particle emission and -1 for local charge conservation. Preliminary NA49 
measurements are plotted in Fig.[| right versus the ratio (N ch ) / (N ch ) tot of the multiplicity 
in the acceptance window and the total multiplicity in the events. Again no significant 
energy dependence is observed and the results are close to the prediction for independent 
particle emission plus global charge conservation = \Jl — (N c h) / (N ch ) tot — 1 (ref. |22[ , 
line in Fig||). 




Figure 10: Fourier coefficients V\ (directed flow) and v% (elliptic flow) of the azimuthal 
distribution of pions in min.bias Pb+Pb collisions versus the rapidity at 40 (left) and 158 
(right) A-GeV beam energy (NA49 preliminary). 
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6 Anisotropic flow 



Anisotropic flow in non-central collisions is sensitive to pressure in the early stage of 
the reaction, which can transform the initial space anisotropy of the reaction zone into 
an azimuthal anisotropy of the momentum distribution of the observed particles. The 
onset of deconfinement might result in a minimum of this effect [f|. Anisotropic flow is 
quantified by the Fourier coefficients v n of the distribution of particle azimuthal angles $ 
with respect to the reaction plane \I/ [23]: 

v n = (cos(n(^i - *„))) = y/2{sin{n ■ • sin(n ■ * n )) 
Corrections for reaction plane resolution, nonuniform azimuthal acceptance and mo- 
mentum conservation were applied. The results for v\ (directed flow) and v 2 (elliptic flow) 
for pions are plotted versus rapidity in Fig.[L0]. The values for v± decrease from 40 (left) 
to 158 (right) A-GeV by a factor of 2, whereas t>2 shows a slight increase. 



7 Conclusion 

The study of central Pb+Pb collisions in NA49 at 40, 80 and 158 A-GeV led to the 
following conclusions: 

• the produced number of pions per participant in Pb+Pb collisions changes from 
suppression with repsect to p+p reactions to enhancement in the SPS energy range 

• the fraction of produced particles containing sors quarks passes through a maxi- 
mum at low SPS energies 

• strangeness production starts to be undersaturated with respect to statistical equi- 
librium at SPS energies at a level consistent with the deconfinement hypothesis 

• no unusual features are found in the evolution of other characteristics of the pro- 
duced hadron system 

NA49 will close the data gap between existing measurements at the AGS and the SPS 
with runs at 20 and 30 A-GeV in 2002. 
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